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SUMMARY 

\  «  - 

review  has-been  made  of  the  basic  shock  tube  along  with  the 
various  modifications  required  to  produce  hypersonic  flow  of  short 
duration.  Modifications  to  the  driver  system  include  multiple 
diaphragms,  area  discontinuities  at  the  diaphragm  station  and 
confcustion  drivers.  The  Influence  of  real  gas  effects,  including 
shock  wave  attenuation,  has-been  noted  for  both  the  production  of 
strong  shocks  (  Mg  >  3  )  and  the  creation  of  hypersonic  flow  in  an 
expansion  nozzle.  Nonreflected,  reflected  and  tailored-interface 
type  shock  tunnels  are  discussed  in-  Section  6  along  with  their 
starting  problems,  real  gas  effects  and  Reynolds  number  and  stag¬ 
nation  temperature  simulation.  Detailed  calculations  have- been 

.  t  • 

omitted  for  simplicity  but  many  figures  have-been  presented  which 
Illustrate  the  operation  of  the  various  shock  tube  configurations. 
Further  details  may  be  obtained  from  the  references  given.  The 
advantages  and  disadvantages  of  the  various  methods  of  producing 
hypersonic  flow  are  summarized. on  the  foilwing  pageG. 
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DRIVER  SY3TENS  FOR  HYPERSONIC  SHOCK  TUNNELS 


Type 


Advantages 


I.  Uniform  Tube, 
Single  Diaphragm 
Station 


D 


ft 

Drfvtr  •* 

—  DWv»m  — — « 

Pi*P> 

2.  Multiple  Diaphragm 
Stations 


P _ JP 


Pa 

* - briv 

•  Ormn- 

1*  Simplest  method 


1«  Overall  pressure  ratio 
reduced  for  a  given  shock 
Mach  number  or  a  gain  in 
Mg  for  a  given  Pj^. 

2*  May  he  combined  vith 
buffer  gas  technique  in 
order  to  use  cold  drivers. 


3.  Area  Discontinuity 
at  the  Diaphragm 

Station,. 

»*> _ 


D 

P* 

M 

PA'P' 

* - M 

vtr - - 

Drhtifi 

Combustion  Driver  • 


1.  Overall  pressure  ratio 
less  than  that  of  (l)  for 
same  conditions. 

2.  May  be  combined,  vith  (2) 
to  give  an  additional  gain 
in  performance  over  (2). 


1.  May  be  used  for  any  of 
the  above  methods. 


I  Ht-$oX 

D 

Hi<  v/»»* 

A 

Hx-Prt  Pa 

Orlw  *n 

Pa  »P, 

2.  Increases  both  pressure 
and  internal  energy  ratio 
across  diaphragn. 

3.  Theoretically,  most  eff¬ 
icient  form  of  driver 


*• 


Disadvantages 

1«  Large  diaphragm 
pressure  ratio  re¬ 
quired  to  produce 
strong  shocks. 


1.  Added  complexity. 
Probably  not  too  prac¬ 
tical  when  more  than 
two  diaphragms  are 
used. 


1.  Diaphragm  pressure 
ratio  required  is 
still  relatively  high 
when  used  vith  cold 
drivers. 

2.  Driver  section 
must  be  reinforced 
due  to  larger  area. 

1.  Added  complexity 
of  igniting  driver  gas. 

2*  May  be  erratic  in 
operation-detonation 
is  a  problem. 

3.  Difficult  to  ach¬ 
ieve  constant  volume 
burning  which  is  most 
efficient  conbustion 
system. 

4.  Attenuation  greater 
than  cold  drivers 

5.  Driver  section  must 
be  reinforced  to  with¬ 
stand  increased  press¬ 
ure. 
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5.  Buffer  Gas 


£ 


D 

* 

.IaI 

A 

-Drivtr- 

-Driven- 

/*>/>**  A 


1.  Enables  cold  hydrogen 

to  be  used  as  driver  vithout 
creating  combustion  at  the 
contact  surface. 

2.  Reduced  overall  diaphragm 
pressure  ratio  for  same  shock 
Mach  number  or  increased  Ms 
for  same  overall  P41. 

3.  Combined  with  (3)  to 
give  performance  equal  to 
(•'0. 


1.  Added  complexity 
of  introducing  a 
third  gas. 

2.  Two  diaphragm 
stations  required. 


4.  Attenuation  problem 
alleviated. 

5.  Extends  range  of  tailored 
interface  method  for  cold 
hydrogen  driver  to  Mg  =  15. 
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HYPERSONIC  TEST  FACILITIES 


Type 

1*  Simple  Shock  Tube 


Advantages 

1.  Simplest  method. 

2.  Simulate  Re  and  temp¬ 
erature  without  any  modif¬ 
ications  . 


Disadvantages 

1.  Maximum  Mach 
number  in  test  region 
is  aproximately 
three. 

2.  A  long  tube  is 
needed  to  obtain  a 
reasonable  test  time. 

3.  Attenuation  pro¬ 
blems. 

4.  Relatively  small 
models. 


2.  Divergent  Expansion 
Nozzle  -  Ndn-reflected 
Method 


1.  Larger  models. 

2.  Increase  in  test  Mach 
number  for  same-  shock  Mach 
number. 

3*  May  be  added  at  exit  of 
type  1  without  modifying  ex¬ 
isting  structure. 


1.  Moderate  Re  sim¬ 
ulation. 

2.  Attenuation  still 
a  problem. 

3.  Tube  must  be  as 
long  as  type  1. 

4.  Running  time  only 
slightly  longer  than 
type  1. 

5.  Starting  problems 
may  require  the  add¬ 
ition  of  a  plenum 
chamber. 

6.  Two  step  nozzle 
may  be  required  for 
hypersonic  Mach 
numbers. 


3»  Reflected  Method 


1.  Large  models. 


2.  Higher  test  Mach  number 
for  same  Ms  than  types  1  and  2. 

3*  Attenuation  problem  only 
affects  production  of  reservoir 
conditions  so  smaller  diameter 
tube  may  be  used. 

4.  Slightly  longer  test  time 
than  types  1  and  2  for  same 
test  Mach  number. 


1.  Low  Re  simulation 

2.  Tube  still  rela¬ 
tively  long. 

3.  For  same  shock 
Mach  number  test  time 
reduced  due  to  refl¬ 
ected  shock  interact¬ 
ing  with  the  contact 
surface. 
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Disadvantages 


Type 

3.  (contd. ) 


5.  Starting  pressure  ratio 
less  than  type  2 


v 

V 


4.  Tailored-Interface  1.  Large  Model. 

Method  -  Unsteady 

Modification  2.  Same  test  Mach  number 

range  as  type  three. 


3.  Test  times  nearly  8  times 
greater  than  type  2. 

k .  Charnel  length  may  be 
reduced. 

5.  Attenuation  less  due  to 
shorter  channel. 


4.  Stagnation  temp¬ 
erature  may  be  high 
enough  to  cause  for¬ 
mation  of  nitric 
oxide  which  may  be 
frozen  by  expansion. 

1.  Tow  Re  as  for 
type  3. 

2.  Running  times 
are  lorig  enough  to 
cause  «rosion  at 
nozzle  throat. 

3.  Radiation  losses 
may  be  important. 

4.  Driver  gas  must 
be  heated  at  higher 
Mach  numbers  in  or¬ 
der  to  achieve 
tailoring  (see  also 
Page  it-; -type  J.) 

5.  Longer  driver 
section,  must  be  long- 
er  than  for  type  3. 

6.  Driven  section 
must  withstand  high¬ 
er  pressure  than 
types  2  and  3. 


5.  Tailored- Interface  1.  Test  times  are  25  times  1.  Nozzle  at  dia- 

Method  -  Steady  that  of  type  2.  phragm  station  must 

Modification  be  changed  for  each 

2.  Other  advantages  as  above.  test  Mach  number. 


2.  Diaphragn  press¬ 
ure  ratio  is  higher 
than  type  4. 

3.  Other  disadvan¬ 
tages  as  above. 


1 — -- 
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1.  INTRODUCTION 

During  the  past  five  years  the  shock  tube,  and  its  various  develop¬ 
ments, have  achieved  prominence  as  a  device  for  the  study  of  hypersonic 
flow.  This  apparatus  has  an  advantage  over  the  wind  tunnel  in  that 
both  Mach  number  and  temperature  can  be  simulated.  There  is  a  disad¬ 
vantage  in  that  total  testing  times  are  very  short,  of  the  order  of 
one  millisecond,  so  that  special  instrumentation  techniques  are  required 
for  the  measurement  of  the  flow  quantities.  However,  methods  have  been 
developed  for  physical  measurements  in  the  short  time  available  and  quite 
acceptable  results  have  "been  obtained. 

It  is  the  purpose  of  this  report  to  outline  some  of  the  more  recent 
developments  of  the  shock  tube  principle.  A3  these  developments  are  only 
modifications,  a  brief  resume  of  the  theory  and  operation  of  the  simple 
shock  tube  will  also  be  presented.  A  more  historical  introduction  is  to 
be  found  in  Refs.  1,  and  2  which  contain  many  references  to  earlier 
research  in  this  field.  A  review  of  more  recent  advances  is  given  by 
Hertzberg  (Ref.  3)  and  Glass  and  Hall  (Ref.  4). 

2.  THE  SIMPLE  SHOCK  TUBE 

The  simple  shock  tube  consists  of  a  straight  tube  of  uniform  cross- 
section  separated  by  a  thin  diaphragm  vhich  divides  the  tube  into  two 
compartments  and  enables  a  pressure  ratio  to  be  created  between  them. 

One  compartment  is  lmovra  as  the  high  pressure  chamber,  or  driver  section, 
the  other,  the  low  pressure  channel,  or  driven  section.  Usually  the 
driver  section  is  pressurised  while  the  driven  section  may  he  either 
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evacuated  or  at  atmospheric  pressure.  The  gases  in  both  sections  are 
normally  in  thermal  equilibrium.  Upon  rupturing  the  diaphragm  the 
wave  system  shown  in  Fig.  1  is  established. 

Theoretically,  upon  removal  of  the  diaphragm,  a  shock  wave  is 
propagated  into  the  low  pressure  section  and  the  high  pressure  gas 
expands  into  the  driven  section  by  means  of  a  rarefaction  wave  cen¬ 
tered  at  the  origin.  This  is  not  quite  true  in  practice  as  the 
rarefaction  wave  is  not  centered  at  the  origin  and  the  shock  wave 
takes  a  finite  time  to  form  from  a  series  of  compression  waves  created 
when  the  diaphragm  is  removed.  However,  these  deviations  from  the 
ideal  theory  are  not  too  important  for  practical  purposes  and  after 
the  shock  has  formed  it  is  considered  to  travel  at  a  constant  velocity. 

The  gas  in  the  driven  section  is  compressed  and  heated  by  the 
shock  wave  while  the  driver  gas  is  expanded  and  cooled  by  the  rarefaction 
wave.  Thus  two  bodies  of  gas  exist  in  the  shock  tube  which  are  brought 
to  the  same  pressure  and  have  the  same  particle  velocity  but  due  to 
different  iermative  processes  their  temperature,  density  and  entropy 
are  different.  These  two  states  are  separated  by  an  interface  or 

f 

contact  surface  which  in  practice  is  more  of  a  region  rather  than  a 
plane  surface.  This  contact  surface,  or  contact  front  as  it  is  called 
by  some  authors,  travels  at  the  particle  velocity.  Therefore,  behind 
the  shock  there  is  a  -region  of  steady  flow  at  a  high  pressure  and 
temperature  while  behind  the  contact  surface  the  flow  is  again  steady 
with  the  nemo  pressure  but  at  a  lower  temperature .  Due  to  the  steady  ■ 
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nature  of  these  states  many  attempts  have  been  made  to  utilize  them 

for  aerodynamic  testing*  The  particle  velocities  are  identical  in 

♦ 

each  state  but  due  to  the  lower  temperature  in  (3)  the  Mach  number 

(M3)  is  higher  than  in  (2)*  However,  state  (2)  is  more  adaptable  to 

hypersonic  testing  due  to  temperature  requirements  and  in  practice 

is  found  to  be  much  more  uniform  than  region  (3)»  ^  ^  ^  The 

testing  time  is  determined  by  the  length  of  the  driver  and  driven 

sections,  the  gas  combination  used  and  the  strength  of  the  initial 

(2) 

shock  wave.  If  the  shock  tube  is  too  short  the  flow  will  be 
temtnated  very  quickly  by  reflected  waves  from  the  dosed  ends. 

Ibr  an  infinitely  long  tube  the  test  time  at  a  station,  X  is  the 
interval  between  the  arrival  of  the  initial  6hock  wave  and  the 
passing  of  the  contact  surface.  Therefore,  for  strong  Bhockp,vVhere 
the  contact  surface  velocity  is  high,  the  test  section  should  be  as 
far  as  possible  from  the  diaphragm  station.  The  total  lengths  of 
the  driver  and  driven  sections  must  be  chosen  so  that  the  reflected 
shock  and  rarefaction  waves  do  not  return  before  the  arrival  of  the 
contact  surface .  The  ideal  lengths  will  be  such  that  the  reflected 
Ghock  and  rarefaction  meet  the  contact  surface  at  the  test  section. 

3.  BASIC  SHOCK  TUBE  EQUATIONS 

The  states  on  each  side  of  the  rarefaction  are  related  by  the 
isentropic  rdationships  while  the  normal  shock  relationships  are 
employed  for  the  transition  across  the  shock  front.  These  two  solutions 
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can  be  matched  at  the  contact  surface  and  the  physical  properties  of 
the  complete  shock  tube  flow  determined.  This  procedure  leads  to  the 
following  identity  which  is  the  basic  equation  for  the  shock  tube. 

'*■- i  ['-ft- ' 

Equation  (3:l)  indicates  that  the  shock  pressure  ratio  (Pg^)  is  a 
function  of  the  diaphragm  pressure  ratio  the  internal  energy 

ratio  (E^)  and  the  specific  heat  ratio  of  the  driver  and  driven 
gases  (  ^  and  ).  The  requirements  for  producing  strong  shock 
waves,  which  are  necessary  in  order  to  achieve  temperatures  suitable 
for  hypersonic  simulation,  may  be  better  examined  by  allowing  the 
diaphragm  pressure  ratio, (Pj^)  to  equal  zero  in  the  above  equation. 
This  gives  the  following: 


•f-  <*  t 

2  G4F/4 


Upon  examination  of  this  equation  it  is  found  that  two  methods  present 
themselves  for  the  production  of  strong  shock  waves  In  a  simple  shock 
tube.  The  diaphragm  pressure  ratio,  (P^)  can  be  made  very  large  or 
the  energy  ratio  across  the  diaphragm^Ej^)  can  be  made  very  small,  or 
both.  The  energy  ratio  may  be  reduced  by  heating  the  driver  and  cool¬ 


ing  the  driven  gases,  but  a  more  cosmon  method  is  to  use  a  light  driver 
gas  and  a  heavy  driven  gas.  The  aerodynamic! st  usually  likes  to  use 


(3*2) 
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air  as  a  test  gaSjWith  possibly  nitrogen  as  a  second  choice, with  helium 
or  hydrogen  drivers.  Figure  2  gives  the  diaphragm  pressure  ratio  re¬ 
quired  for  a  range  of  shock  pressure  ratios  for  the  combination  Air/Air, 
He/Air,  and  l^/Air  for  the  simple  shock  tube  case  while  Fig.  3  presents 
the  variation  of  shock  wave  Mach  number, (Mg)  and  particle  velocity, (Ug^) 
with  diaphragm  pressure  ratio ,(Pj^)  for  the  three  gas  combinations. 

It  can  be  seen  from  these  two  figures  that  a  great  advantage  is 
gained  by  the  use  of  helium  or  hydrogen  as  a  driver  gas.  The  internal 
energy  ratio  can  also  be  reduced  by  heating  the  driver  gas  either  by 
electrical  heating  or  by  the  ignition  of  combustible  mixtures  which 
give  high  values  of  T4.  The  combustible  mixture  is  usually  a  stoichio¬ 
metric  mixture  of  hydrogen  and  oxygen  diluted  with  excess  hydrogen  or 
helium.  ^  This  method  has  seen  widespread  use  for  the  production  of 
strong  shock  waves  in  shock  tubes.  Another  method  used  in  gun  tunnel 
versions  of  the  shock  tube  is  to  compress  the  driver  gas  by  a  piston 
driven  by  an  explosive  process,  usually  the  ignition  of  a  powder  charge* 
This  compression  raises  the  pressure  and  temperature  to  a  predetermined 
level  governed  by  the  strength  of  the  diaphragm  which  then  bursts 
initiating  the  flow. 

Having  established  Eq.  (3:l)  all  the  other  flow  quantities  can  be 
determined.  The  most  useful  equations  are  listed  below  in  nondimension- 
al  notation.  ^ 
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1.  Density  ratios 


(3:3) 


r. 


21 


I 


«>  *  n, 

2.  Speed  of  sound  and  temperature  ratio 

A34  —  T34  -  ~  \pH  /J^| 


A%i  ss  %!  ss 


L  i  Rti  J 

3.  Velocity  of  the  shock  wave 


Pz\  (°*l  *  ft/) 


Ms  ^  0  +  >?/)J 


J. 
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4.  Particle  velocity  or  contact  surface  velocity 

Mu  _  Pn  ~i 

X  [(/,(*,  P„+iJf 

M»«  -  —L.  [/  -  (P^P,,)*]  =  AnU, 


5.  Local  Mach  numbers 


VI 
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)h± 
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(3:5) 
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(3:10) 
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6.  Speed  of  the  head  and  tail  of  the  rarefaction  wave 


^34  A  34 


(3:12) 


-/ 

The  head  of  the  rai ei' action  wave  travels  at  the  velocity  of  sound  of 
the  driver  gas. 

It  is  of  interest  to  examine  the  limits  of  these  equations  for  an 
infinite  value  of  the  diaphragm  pressure  ratio.  Table  I  is  a  compari¬ 
son  of  three  different  driver  gases,  air,  helium,  and  hydrogen  in  com¬ 
bination  with  air  as  a  driver  gas.  Table  I  indicates  more  than  a  ten¬ 
fold  increase  in  shock  strength  when  hydrogen,  rather  than  air,  is  used 
as  a  driver  gas.  This  same  increase  is  obtained  in  the  temperature 
achieved  in  the  region  behind  the  shock,  however,  the  Mach  number  in 
the  steady  state  region  (2)  is  limited  by  this  temperature  increase, 
and  even  for  the  ultimate  case,  where  P-^  and  are  both  zero,  the 
limiting  Mach  number  (M2)  is  only  I.89.  Thus  Table  I  illustrates  that 
although  high  temperatures  may  be  achieved  this  Mach  number  restriction 
prevents  true  hypersonic  simulation  in  a  simple  shock  tube.  The  results 
shown  in  Table  I  are  based  on  perfect  gas  theory  (i.e.  Y  «t  1,4)  which 
does  not  hold  for  the  extreme  shock  pressure  ratios  calculated. 


4.  MODIFICATION  OF  THE  SIMPLE  SHOCK  TUBE  FOR  THE  PRODUCTION  OF  STRONG  SHOCK  WAVES 
4.1  INTRODUCTION 

It  is  possible  to  improve  the  performance  of  the  simple  shock  tube 
by  suitable  modifications.  These  modifications  include  the  installation 
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or  additional  diaphragm  stations  and  driver  sections  whose  cross-seetional 
area  is  greater  than  the  drivdn  section. The  former  method  is 
a  device  which  sacrifices  pressure  ratio  in  favor  of  an  increased  temp¬ 
erature  ratio.  The  latter  creates  an  additional  expansion  at  the  diaphragm 
station  which  converts  thermal  energy  into  kinetic  energy.  Therefore, 
in  both  cases,  it  is  possible  to  produce  the  same  shock  strength  as  a 
uniform  tube  but  with  a  smaller  diaphragm  pressure  ratio,  or  conversely, 
to  produce  a  stronger  shock  for  the  same  diaphragm  pressure  ratio. 

Another  method,  which  increases  the  temperature  ratio  across  the  dia¬ 
phragm,  is  to  use  a  combustible  mixture  as  a  driver  gas.  Rather  than 
use  combustible  mixtures  it  is  sometimes  more  satisfactory  to  use  either 
cold  helium  or  hydrogen  as  a  driver,  as  mentioned  in  Section  (3).  A 
further  modification  of  this  technique  has  been  to  introduce  an  inert 
buffer  gas  between  the  hydrogen  driver  and  the  air  test  gas  with  an 

(12)  (13) 

area  discontinuity  at  the  buffer-air  diaphragm  station.  These 

methods  will  now  be  examined  in  more  detail. 

4.2  MULTIPLE  DIAPHRAGM  SHOCK  TUBES 

Figure  4  is  a  sketch fof  a  shock  tube  of  constant  area  utilizing 

two  diaphragms.  D^,  is  the  primary  diaphragm  which  is  ruptured  first 

producing  the  shock  5g.  This  shock  wave  is  reflected  at  diaphragm 

Dg  bringing  the  gas  in  state  5  to  rest  and  raising  its  temperature 

end  pressure.  After  a  short  delay,  diaphragm  D2  is  removed  and  a 

new  wave  system  is  produced  resulting  in  the  shock  wave  S^  which  is  of 

greater  strength  than  the  shock,  Sg.  Comparing  Fig.  1  and  Fig.  4  for  * 

the  same  overall  temperature  and  pressure,  and  with  T-j_=T£=Tg,  U  Vg 

'•J 
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will  give  an  idea  of  the  Increase  in  shock  strength  for  the  double 
diaphragm  case.  As  the  gas  is  brought  to  rest  ahead  of  Dg  the  pressure 
and  temperature  ratios  of  interest  in  both  cases  are  and  Tj^. 
However,  Pj^<  P^  while  >Tjj  which  gives  a  shock  of  greater 
strength  than  the  system  shown  in  Fig.  1.  As  an  example  of  the  gain 
to  be  ejected  by  the  addition  of  diaphragm  Dg  consider  the  case  where 
the  overall  pressure  ratios,  P^  (Fig.  l)  and  P81  (Fig.  4)  are  both 
equal  to  10^.  Then  from  the  Air/Air  curve  of  Fig.  3>  Mg=3*l4.  Fbr  the 
case  shown  in  Fig.  4  let  Pgj_=10  which  means  pQg=102.  The  initial 
shock  wave  Mach  number,  Mjg  will  be  found  using  Fig.  3  and  is  2.37. 

The  final  overall  pressure  ratio  will  be  P^  a  p^g  x  Pg-j_.  The  value 
of  Pj;g  is  obtained  from  interpolation  of  Table  III  in  Ref.  (9) 

(p46  -  25.85)  while  T^g  «s  3«236.  Having  determined  these  new  initial 
conditions,  the  final  shock  Mach  number^  (*3.77)  con  be  obtained  by 
interpolation  from  Table  I  Ref.  (9).  It  is  seen  that  a  gain  of  20$ 


in  shock  Mach  number  is  obtained  by  the  addition  of  the  diaphragm,!^. 

The  intermediate  pressure,  Fg  can  be  adjusted  so  as  to  give  a  maximum 
Ghock  Mach  number  for  a  given  overall  pressure  rctio,  p8r 

The  above  method  may  be  somewhat  simplified  by  making  Dg  weak  enough 


that  the  initial  shock  vave,Mg£  will  rupture  it  on  contact.  This  produces 
on  unsteady  expansion  from  state  6  to  state  1  resulting  in  a  less  efficient 


process  than  the  reflected  shock  type  described  above,  but  reducing  the 


complexity  of  the  system  caused  by  the  introduction  of  a  mechanical  or 


electrical  delay.  The  vavo  diagram  for  this  type  of  shock  tube  is  shown 
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in  Fig.  5»  The  onalysis  for  this  case  is  carried  out  in  Rei’. 

Carrjing  out  the  some  calculation  for  the  same  conditions  as  used  for 
the  example  of  the  reflected  shock  type  of  double  diaphragm  shock  tube 
gives  a  final  Mach  number, of  3.73  which  is  only  a  decrease  of  app¬ 
roximately  one- percent.  This  variation  would  be  larger  at  higher  Mach 
numbers  but  the  percentage  decrease  would  probably  never  be  large  enough 
to  warrant  adopting  the  delay  type  installation  for  a  straight  shock  tube. 

Ref.  9  indicates  that  for  the  maximum  gain  in  Mg  the  pressure  must 
be  approximately  the  geometric  mean  of  the  pressures  on  either  side. 

Then  for  a  multiple  diaphragm  shock  tube  the  overall  pressure  ratio  is 
given  by  PQ  =  Pj_n,  where  is  the  individual  pressure  ratio  across 
each  diaphragm  and  n  is  the  number  of  diaphragms.  It  can  be  shown 
(Ref.  9)  that  there  is  a  maximum  shock  Mach  number  attainable  irres¬ 
pective  of  the  number  of  diaphragms  used  and  depending  only  on  the 
overall  pressure  ratio,  across  the  shock  tube.  This  relationship  is 
Mg  Max.  =  P0  3/1^  and  Table  II  gives  the  values  of  Mg  Max.  attainable 
over  a  range  of  reasonable  overall  pressure  ratios  for  the  conditions 
of  equal  and  T  throughout.  It  is  also  shown  in  Ref.  9  that  a  multiple 
diaphragm  system,  where  all  the  diaphragms  except  the  first  have  zero 
pressure  difference  across  them,  is  another  method  for  creating  strong 
shock  waves.  If  a  chock  wave  of  Mach  number  greater  than  2.67  is 
generated  by  rupturing  the  first  diaphragm  then,  theoretically,  the 
maximum  obtainable  shock  speed  increases  without  limit.  In  this  system 
the  shock  wave  reflecting  off  each  successive  diaphragm  will  increase 
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the  pressure  and  temperature  ahead  of  that  diaphragm  producing  a  succ¬ 
essively  stronger  chock. 

4.3  AREA  DISCOHTIHUITIES 

As  outlined  in  Section  2  the  basic  shock  tube  comprises  a  driver 
and  driven  section  of  uniform  cross-section  separated  by  a  diaphragm. 
However,  many  authors  have  shown  that  there  is  a  gain  in  the  final  shock 
strength  by  introducing  an  area  discontinuity  at  the  diaphragm  station. 
This  arrangement  may  also  be  combined  with  the  double  diaphragm  method 
to  give  a  further  increase  in  the  final  shock  strength.  This  latter 
method  is  illustrated  in  Fig.  6.  The  area  change  may  also  be  accom¬ 
plished  by  a  transition  section  as  shown  in  Fig.  7. 

In  Fig.  6  if  diaphragm  Dg  is  eliminated  a  single  shock  system  is 
set  up  at  Djlj  with  being  replaced  by  a  rarefaction  wave  as  in  the 
case  of  a  uniform  tube.  If  the  shock  S-^,  is  such  that  1  then  the 
wave  system  will  be  as  shown  in  Fig.  1,  the  flow  being  accelerated 
subsonically  to  However,  if  M3  > 1  the  gas  in  (6)  is  accelerated  to 
sonic  velocity  at  the  diaphragm  station  and  is  then  further  expanded  by 
a  rarefaction  wave  to  the  supersonic  Mach  number,  M^.  The  addition  of 
the  diaphragm,  Dg  combines  the  multiple  diaphragm  and  variable  geometry 
type  using  a  convergent  transition  section  which  may  or  may  not  contain 
a  convergent  -  divergent  nozzle.  However,  the  only  advantage  here  is  a 
smOoVik  r  flow  from  the  driver  to  the  driven  section  as  the  flow  processes 


arc  ip'-’ 


’hither  the  area  change  is  gradual  or  discontinuous. 
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As  in  the  previous  section  a  comp orison  will  he  made  between  the 
uniform  and  the  variable  geometry  shock  tubes  by  determination  of  the 
final  shock  Mach  number  achieved  for  identical  overall  pressure  ratios. 
An  analytical  expression  for  such  a  gain  factor  is  not  possible  but  by 
use  of  tables  and  graphs  it  is  possible  to  determine  these  values  of 
•  We  will  only  consider  the  case  of  which  corresponds  to  the 
range  of  shock  strength  required  for  hypersonic  testing.  Consider  first 
the  case  where  the  diaphragm,  Dg  is  not  present  and  that  A^  =A^  is  the 
minimum  section  such  as  shown  in  Fig.  7*  Reference  (10)  develops  the 
following  expression  for  the  »  -erall  pressure  ratio,  P^. 


a- 


ft 

9 


l  +  hi!  M, 
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where  9  is  an  equivalence  factor  given  by 
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As  is  the  minimum  section,  will  be  equal  to  unity  but  Mg&  must 
be  determined  in  order  to  calculate  9  »  Knowing  the  area  ratio,  Ajt/A^ 
then  Mo  is  determined  from  the  following  relationship. 


(4:3:1) 


(4:3:2) 
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It  is  seen  from  Eq[S.  (4:3:2)  and  (4:3:3)  that  for  a  given  area 
ratio  5  is  a  constont  when  >1*  Therefore  knowing  the  area  ratio, 
g  can  he  determined  and  a  plot  of  vs.  Pg^  can  he  obtained.  A 
comparison  of  the  shock  Mach  numbers  attainable  for  area  ratios  of  1, 
and  *6  is  presented  in  Fig.  8e  Table  3  is  token  from  Ref.  (10)  and 
gives  the  maximum  shock  Mach  numbers  when  P^i"*00  for  various  gas  com¬ 
binations,  four  area  ratios  and  two  temperature  ratio.  Both  Fig.  8 
and  Table  3  illustrate  that  the  gain  in  shock  strength  is  not  too 
great  even  when  the  area  ratio  is  infinite.  Once  again,  however,  the 
gain  achieved  by  increasing  the  sound  speed  in  the  driver  gas  is  in¬ 
dicated  by  Table  3* 

As  an  example  of  the  gain  to  be  expected,  through  the  use  of  an 
area  change  at  the  diaphragm,  consider  the  example  used  previously. 

In  this  Instance  the  area  ratio  will  be  3:1,  the  diaphragm  pressure, 
p4l  *  lo3  for  an  Air/Air  combination  with  a  constant  initial  temperature 
ratio,  Tj^,  which  gives  a  final  shock  Mach  number  of  3*35*  Combining 
the  multiple  diaphragm  method  (two  diaphragms)  with  the  same  area 
discontinuity  and  total  pressure  ratio  and  with  Pg£  (Fig.  6)  =  10p 
gives  Mg  za  3.9,  This  Is  on  increase  of  16  per  cent  due  to  the 
addition  of  Dg  to  the  system,  and  an  increase  of  rrarly  25  per  cent 
over  a  simple  shock  tube  with  the  some  Pj^. 

The  above  example  is  for  an  air  driver  but  we  have  seen  from  Fig.  3, 
and  Tables  1  and  3,  that  a  light  gas  such  as  hydrogen  or  helium  is  a  more 
efficient  driver.  However,  when  hydrogen  is  used  to  produce  shocks 
cr.  '-nation  occurs  nt  the  dviver-driven  interface  which  sets  up  severe 
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flow  distui  l'-'nces.  Helium  could  "be  substituted  as  a  driver  but  it  is 
not  us  efficient.  A  method  to  eliminate  this  combustion  problem  is 
the  use  of  a  double-diaphragm  shock  tube  with  an  inert  buffer  gas 
between  the  driving  hydrogen  and  the  driven  air.  There  is  also  an 
area  change  incorporated  at  the  diaphragm  separating  the  buffer  gas 
and  the  air.  Figure  9  shows  the  effect  of  varying  the  pressure  and 
molecular  weight  in  various  buffer  gases  for  an  overall  pressure  ratio 
of  lo\  It  is  seen  from  the  figure  that  a  shock  Mach  number  of  approx¬ 
imately  15  can  be  obtained  by  the  use  of  any  of  the  inert  gases  given 
on  the  figure,  however,  the  lighter  gases  such  as  helium  require  a 
large  pressure  ratio  across  the  downstream  diaphragm.  Consider  as 
an  example  a  shock  tube  of  this  type  with  argon  as  the  buffer  gas.  In 

order  to  produce  a  shock  Mach  number  of  13  a  simple  shock  tube  using 

k 

hydrogen  as  a  driver  would  require  a  Pj(1  of  7  X  10  .  If  a  large  area 
discontinuity  is  introduced  the  diaphragm  pressure  ratio  is  reduced  to 

Ij. 

1.8  X  10  .  When  an  argon  buffer  is  introduced  between  the  driver 
hydrogen  and  driven  gas  an  overall  pressure  ratio  of  only  3*5  X  10^ 
is  required  which  is  a  reduction  by  a  factor  of  20.  This  method  seems 
to  offer  very  good  possibilities  of  producing  strong  shock  waves  using 
cold  driver  gases. 

4.4  C0M3USTI0N  DRIVERS 

A  common  method  of  producing  strong  shock  waves  is  by  using  a 
combustible  mixture  which  is  ignited  to  give  high  pressures  and  temp¬ 
eratures  in  the  driver.  A  typical  composition  for  a  combustible  driver 
is  13.33  per  cent  hydrogen  6.67  per  cent  oxygen  and  80  per  cent  helium 
as  a  diluent.  The  ignition  of  the  hydrogen  and  oxygen  raises  the 


CONVAIR  -  San  Diego 


Report  No.  ZR  658-050 


temperature  and  pressure  of  the  helium  which  gives  a  large  increase  in 
the  energy  and  pressure  ratios  Across  the  diaphragm.  Excess  hydrogen 
has  also  "been  used  as  a  diluent.  Using  this  type  of  driver  the  diaphragn 
is  usually  hurst  by  the  combustion  of  the  mixture.  Combustion  may  be  produced 
in  three  ways.  A  spark  pD.ug,  or  plugs,  may  be  inserted  in  the  driver 
section  -which  is  fed  by  a  high  voltage  coil,  producing  the  spark  which 
causes  ignition.  A  heating  wire  may  be  inserted  the  length  of  the  high 
pressure  chamber  and  a  large  current  passed  through  it  which  causes  a 
constant  volume  burning.  The  mixture  may  be  allowed  to  ignite  itself 
by  introducing  the  combustible  gases  and  then  adding  the  diluent  until 
the  diaphragm  ruptures,  whereupon  the  hydrogen-oxygen  mixture  will  ignite. 

It  has  been  found  that  the  constant  volume  process  is  the  most  efficient. 

There  is  also  a  practical  limit  to  which  the  driver  gas  can  be  heated  due 
to  molecular  dissociation  at  the  higher’  temperatures.  Figure  10  compares 
the  theoretical  shock  Mach  numbers  obtained  over  a  range  of  diaphragm 
pressure  ratios  for  helium,  hydrogen  and  combustion  drivers.  Curves  are 
presented  for  both  the  simple  shock  tube,  (A^/Ap=i)  and  one  where  A^/Ap=6.25. 
DEVIATIONS  FROM  IDEAL  FIUID  THEORY 
5.1  REAL  GAS  EFFECTS 

In  the  previous  chapters  the  discussions  and  calculations  have  been 
carried  out  assuming  that  the  specific  heat  ratio,  $  remained  constant. 
However,  for  shock  Mach  numbers  greater  than  three  the  temperatures  pro¬ 
duced  are  high  enough  to  invalidate  this  assumption  and  the  variations 
from  an  ideal  gas  must  be  considered.  Above  Ms  =  3  fairly  high 

temperatures  will  be  achieved  in  state  (2)  and  very  low  temper?  bures  in 
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state  (3)*  For  an  Mg  of  4.0, with  an  Air/Mr  combination,  the  temperature 
Tg  will  he  40°K  which  is  below  the  boiling  point  of  nitrogen  and  oxygen. 
In  state  (2)  the  temperature,  Tg  will  be  over  1200°K  on  the  basis  of 
the  ideal  theory.  At  these  temperatures  the  specific  heat  ratio,  can 
no  longer  be  taken  as  1.4  and  above  Mj  =  3  the  gas  can  no  longer  be 
considered  an  ideal  fluid  and  the  real  gas  effects  must  be  considered. 

We  are  concerned  primarily  with  state  (2)  behind  the  shock,  ns 
this  will  be  ouj*  test  gas,  so  our  discussion  will  be  limited  to  this 
region.  The  internal  energy  content  of  a  diatomic  gas  is  made  up  as 
follows: 

1)  The  kinetic  energy  of  translation  of  a  molecule 

2)  The  energy  of  molecular  rotation  . 

3)  The  energy  of  molecular  vibration 

4)  The  energy  of  molecular  dissociation  into  atomic  groups 

5)  The  energy  of  electronic  excitation 

6)  The  energy  of  ionization 

As  the  temperature  is  increased  the  higher  rotational  and  vibrational 
levels  are  excited  absorbirg  heat  in  the  process.  This  causes  an 
increase  in  the  heat  capacity  which  in  turn  increases  the  specific 
heats,  Cp  and  Cy.  However,  as  they  both  increase  (about  the  same 
amount  then  the  ratio,  K  will  be  smaller.  Added  to  this  effect  the 
specific  heats,  and  the  molecular  weight,  become  pressure  dependent. 

Above  20C0  -  3'X)0°K  dissociation  occurs  giving  nitrogen,  oxygen  end 
their  atomic  states.  These  temperatures  arc  also  high  enough  to  cause 
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chemical  reactions  such  as  the  formation  of  nitric  oxides.  At  still 
higher  temperatures  (  6000°K)  electronic  excitation  and  ionization 

occur  and  the  air  becomes  almost  completely  dissociated  giving  atomic 
oxygen  and  nitrogen  as  before,  small  percentages  of  the  various  oxides 
of  nitrogen,  electrons,  negative  Og  ions  and  positive  ions  of  nitric 
oxide.  These  conditions  would  be  the  same  as  those  produced  at  the 
stagnation  point  of  an  insulated  blunt  body  traveling  at  Mach  number 
of  1 6  at  an  altitude  of  60,000  feet. 

The  transition  through  a  strong  shock  wave  can  be  qualitatively 
described  as  follows.  At  the  shock  front  the  temperature  jumps  to  a 
high  value  as  active  degrees  of  freedom,  translation  and  rotation  are 
excited  almost  immediately  in  the  distance  of  a  few  mean  free  paths. 
After  this  initial  jump  thermal  equilibrium  will  be  reached  in  an. 
exponential  manner  as  energy  is  transferred  to  the  remaining  degrees 
of  freedom  during  the  relaxation  period.  Therefore,  besides  the 
variation  in  the  flow  properties,  some  consideration  must  be  given  to 
the  relaxation  time  to  insure  that  equilibrium  conditions  exist  for 
the  tost  to  be  performed. 

Figure  11  indicates  the  variation  of  the  temperature  Tg  with 
shock  Mach  number,  Mg  for  two  values  of  the  initial  pressure.  It 
is  seen  that  the  real  gas  temperatures  are  lower  at  the  higher  Mach 
numbers  os  expected 'from  the  precceding  discussion.  Figure  12  compares 
the  density  ratio  for  a  real  and  an  ideal  gas  over  the  same  range  of 
l!s  as  used  in  Fig.  11.  Here  the  real  gas  density  is  higher  than  the 
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ideal  gas  case  which  approaches  a  limit  of  6  while  at  Mg  =  12  the  real 
gas  density  ratio  continues  to  increase.  Due  to  the  increase  of  density 
and  decrease  of  temperature  the  pressure  is  found  to  he  very  near  the 
same  for  the  real  gas  as  for  the  ideal  gas.  Therefore,  Equation  3*1 
will  give  a  fairly  good  approximation  for  real  gas  calculations. 

The  decrease  in  temperature  behind  the  normal  shock  in  a  real  gas 
is  more  pronounced  than  the  change  in  specific  heat  ratio  or  the  decrease 
in  molecular  weight.  Therefore,  the  speed  of  sound  in  a  real  gas  is  lower 
for  the  same  shock  strength  than  for  the  ideal  case  giving  a  higher 
attainable  Mach  number,  Mg  behind  the  shock.  This  is  illustrated  in  Fig. 

13.  The  maximum  value  of  Mg  far  the  ideal  case  is  I.89  but  Fig.  13  indi¬ 
cates  that  an  Mg >3  can  be  obtained  for  a  real  gas.  As  the  Mach  number  of  the 
flow  in  region  (2)  is  increased  while  the  temperature  and  speed  of  sound  are 
decreased  then  the  particle  veiod.  ty  must  be  changed  very  little  byihe 
real  gas  effects.  This  is  found  to  be  the  case  with  the  real  values  of 
Ug  being  only  slightly  greater  than  the  ideal  value. 

The  real  gas  curves  shown  in  Figs.  11,  12  and  13  were  computed  using 
the  charts  and  tables  contained  in  references  18  and  22.  The  crossover 
of  the  two  real  gas  curves  in  Fig.  13  is  due  to  the  identical  behaviour 
of  the  speed  of  sound  ratio  with  increasing  temperature .  The  ideal  gas 
curves  were  computed  from.  Equations  3  ;  3*6  and  3°H*  References  18  and  22 

are  two  examples  of  the  many  excel. lent  sets  of  tables  and  graphs  available 
for  the  determination  of  the  flow  properties  behind  strong  shock  waves. 0-6-2$) 
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Table  III  indicated  that  there  was  no  appreciable  gain  in  perfor¬ 
mance  by  introducing  an  area  change  at  the  diaphragm  station.  However, 
these  calculations  were  based  on  the  driver  gas  being  ideal.  It  has 
been  shovn  that  if  the  driver  gas  density  is  inci eased  to  a  very  high 

value  then  there  is  a  considerable  gain  in  performance  by  using  a  large 

(2*0 

driver-driven  area  ratio.  This  difference  is  due  to  the  behaviour 

of  the  attractive  and  repulsive  interraolecular  forces*.  At  normal  den¬ 
sities  these  forces  are  negligible  and  the  performance  of  a  uniform  and 

non-uniform  tube  are  nearly  identical.  However,  if  the  driver  density 

/ 

is  increased  to  a  very  high  value  there  is  a  considerable  advantage  in 
increasing  the  chamber  area  as  the  repulsive  forces  then  predominate. 

It  was  found  in  Ref.  2h  that  an  area  ratio  of  five  was  nearly  equivalent 
to  an  infinite  area  ratio,  while  a  chamber  twice  the  area  of  the  channel 
would  give  half  the  performance  promised  by  the  infinite  value  for  the 
case  of  nitrogen  and  hydrogen  drivers  at  pressures  of  6000  and  2200 
atmospheres  respectively.  It  was  found  that  at  intermediate  densities 
the  performance  relative  to  an  ideal  gas  would  be  almost  the  same  for 
the  non-uniform  tube  and  considerably  reduced  for  the  uniform  tube. 

The  test  time  for  identical length  chambers  is  reduced  when  using  high 
density  drivers  due  to  the  increase  in  sound  speed,  a^  which  means  that 
the  rarefaction  waves  reflected  from  the  end  at  the  high  pressure  section 
will  travel  at  a  higher  velocity.  Increasing  the  chamber  length  would 
recover  this  loos  in  test  time. 
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1? .  2  SHOCK  WAVE  ATTEl'HJATION 


The  ideal  shock  tube  theory  predicts  that  when  the  diaphragm  breaks 
a  flow  will  be  set  up  as  shown  in  Fig.  1.  According  to  this  theory  a 
plane  shock  wave  is  produced  which  travels  down  the  tube  at  a  constant 
velocity  followed  by  two  steady  states  separated  by  a  contact  surface. 

In  practice  however,  this  model  is  not  realized.  ^  Instead  the  shock 
wave  velocity  is  found  to  decrease  with  distance  from  the  diaphragm. 

In  fact  the  shock  does  not  reach  its  calculated  speed  until  a  short 
distance  from  the  diaphragm  due  to  the  formation  process.  The  shock 
only  attains  the  ideal  value  for  weak  shocks,  for  as  the  diaphragm  press¬ 
ure  is  increased  the  shock  never  reaches  the  theoretically  predicted 
speed,  rising  to  a  maximum  some  distance  from  the  diaphragm  and  then 
attenuating  as  it  progresses  down  the  tube.  This  effect  is  also  more 
pronounced  as  the  shock  tube  diameter  is  decreased. 

Due  to  this  shock  attenuation ^he  flow  quantities  in  states  (2) 
and  (3)  are  not  steady  but  vary  with  time  as  the  shock  progresses  down 
the  tube  thereby  making  testing  in  these  regions  very  difficult.  The 
problem  of  shock  attenuation  has  received  a  great  deal  of  attention  in 
the  last  five  years.  ^  Many  of  these  investigators  have  focussed 

their  attention  on  the  viscous  effects  which,  of  course,  are  not  allowed 
in  the  ideal  theory.  There  have  been  studies  made  of  the  boundary 


layer  both  theoretical  and  experimental  to  determine  the  causes  of  the 
(?-9>  30-32) 

attenuation.  A  satisfactory  explanation  has  not  been  arrived 

at  although  the  theories  of  Trimpi  and  Cohen  and  Kirels  and  Dram  give 
good  agreement  at  lower  Mach  numbers.  There  is  probably  come  influence 
a th  :  dir*  d>;\- ;■  1  static  1  r.r.-rjcinV.y  at  the  higher  shock  strengths  where 
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metal  diaphragms  are  used  and  take  a  finite  time  to  open.  This  effect 
may  he  instrumental  in  reducing  the  velocity  which  muGt  match  Ug  at 
the  contact  surface.  The  complete  flow  can  he  considered  to  he  gover- 
ened  hy  conditions  at  the  contact  surface  therefore  a  reduction  in  the 
velocity  Ug  will  produce  a  reduction  in  the  shock  speed. 

The  detailed  processes  of  attenuation  are  not  too  important  in  the 

operation  of  a  shock  tunnel  hut  it  is  deBirable  to  have  as  little  attenuation 

as  possible  so  that  the  flow  conditions  will  he  fairly  steady  during 

the  test  period.  Considerable  work  has  been  carried  out  to  determine 

(28  ) 

what  initisl  conditions  will  give  the  least  attenuation.  '  in 

addition  to  the  causes  given  above  it  has  been  generally  found  that  a 
greater  attenuation  is  experienced  as  the  sound,  speed  of  the  driver  gas 
is  increased.  Figure  14  shows  typical  attenuation  measurements  for 
several  driver  gas  compositions.  The  results  for  helium  give  the  least 
attenuation  but  this  driver  gives  the  lowest  initial  shock  strength. 

The  combustion  driver,  which  is  the  most  efficient  type,  produces  the 
highest  initial  shock  strength  but  at  L/D  of  150  the  shock  has  attenuated 
to  a  lower  speed  than  the  helium  driven  shock.  It  is  necessary  in  both 
straight  shock  tubes  and  non-reflected  type  shock  tunnels  to  have  as 
long  a  driver  section  as  possible  in  order  to  have  reasonable  test 
times.  The  attenuation  could  be  alleviated  by  increasing  the  diameter 
thus  reducing  the  L/D  ratio  but  there  is  a  limit  to  this  due  to  structural 
and  handling  reasons. 

Experimental  eveidence  has  shown  that  the  attenuation  is  inversely 
proportional  to  the  shock  tube  Reynolds  number^a^/y^).  For  a 
given  shock  tube  using  air  ns  the  test  gas  the  Reynolds  number  can 
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only  be  adjusted  through  the  kinematic  viscosity  or  more  precisely 
by  adjusting  the  initial  pressure,  p-^.  From  the  attenuation  evand- 

0 

point,  it  is  desirable  to  have  p^  as  large  as  possible.  However,  if 
cold  drivers  are  used,  which  seems  desirable  from  a  study  of  Figure  14, 
a  low  initial  pressure  will  be  required  to  produce  the  diaphragm 
pressure  ratio  to  give  the  same  shock  strength.  These  requirements 
are  at  cross  purposes  and,  as  usual,  a  suitable  compromise  must  be 
found  in  order  to  reduce  the  attenuation  effects, 
b.  METHODS  OF  GENERATING  HYPERSONIC  FLOW 
6.1  EXPANSION  NOZZLES 

It  has  been  shown  in  the  previous  section  that,  although  the 
proper  temperature  may  be  achieved  in  a  straight  shock  tube,  the 
maximum  Mach  number  is  too  low  for  hypersonic  simulation.  In  ..order 
to  overcome  this  difficulty  the  hypersonic  shock  tunnel  was  developed. 
By  placing  and  expansion  nozzle  at  the  downstream  end  of  the  conventional 
shock  tube  the  flow  behind  the  shock  can  be  expanded  to  a -higher  Mach 

(34-37) 

number.  Therefore,  by  a  proper  choice  of  shock  strength  and  expan¬ 

sion  ratio  it  is  possible  to  simulate  both  temperature  and  Mach  number. 
Figure  15  presents  the  shock  tube  area  ratio  required  over  a  range  of 
test  section  Mach  numbers. 

The  simplest  type  of  installation  is  a  single  stage  expansion 

nozzle  placed  at  the  exit  of  the  driven  section  of  the  tube,  the 

final.  Mach  number  being  determined  by  the  area  ratio.  This  system 

is  known  as  the  nonreflocted  method.  In  Ref.  (35)  a  set  of  mahogany 

liners  were  installed  in  the  shock  tube  with  the  nozzle  located  at  the 

0 
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existing  test  section.  Using  thiG  simple  scheme  a  study  was  made  of 
the  flow  in  the  nozzle.  Figure  l6  illustrates  the  wave  system  set  up 
in  a  shock  tube  of  this  type.  As  the  initial  shock  wave,  enters  the 
nozzle  it  weakens  locally  and  develops  curvature  which  propagates  from 
the  wall  to  the  centerline  of  the  tube  as  the  shock  moves  into  the 
test  section.  As  the  shock  weakens  the  velocity  of  the  particles 
behind  decree  'es,  however,  the  particles  originally  set  in  motion  by  the 
passage  of  the  shock  down  the  uniform  tube,  (Ag)  accelerate  as  they  pass 
through  the  expansion  waves  at  the  mouth  of  the  nozzle  and  increase 
in  velocity.  Therefore,  there  is  another  shock  created  which  separates 
the  two  flow  regions  2  and  4,  (Fig.  16).  Due  to  the  entropy  being 
different  behind  each  shock  there  is  also  a  contact  surface  formed 
which  separates  the  regions  2  and  3 t  (Fig.  16).  The  wave  Sg  attempts 
to  propagate  upstream  but  as  it  is  moving  against  a  supersonic  flow  it 
Is  swept  downstream  through  the  nozzle.  As  S.  moves  through  the  nozzle 
it  weakens  which,  in  turn,  causes  S2  to  become  stronger.  If  83  grows 
sufficiently  strong  it  will  not  be  swept  out  of  the  nozzle  but  will 
remain  standing  in  the  nozzle  preventing  the  tunnel  from  starting. 

Shock  tunnels  have  been  built  which  use  double  and  triple  expansion 
nozzles  for  the  generation  of  hypersonic  flow.  This  system  enables  the 
test  section  area  to  be  reduced  over  that  required  for  the  system  shown 
in  Fig.  16.  In  order  to  obtain  reasonable  testing  times  at  hypersonic 
Mach  numbers  it  is  necessary  to  have  a  fairly  long  shock  tube.  However, 
with  long  tubes  the  boundary  layer  built  up  on  the  walls  may  be  quite 
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thick  causing  severe  attenuation  problems.  A  method  of  circumventing 
this  is  to  increase  the  diameter  but  this  in  turn  means  increasing  the 
test  section  area  of  the  driven  section.  By  introducing  a  two  step 
nozzle,  which  givea  two  stage  expansion  and  also  drains  off  the  boundary 
layer  created  behind  the  initial  shock  wave,  the  test  section  area  can 
be  decreased.  This  type  of  shock  tunnel  is  illustrated  in  Fig.  17* 

6.2  STARTING  PROBLEMS 

The  flow  duration  is  the  time  between  the  departure  of  the  start¬ 
ing  wave  system  S-^  Cg>  Sg,  Fig.  16,  from  the  test  section  and  the 
arrival  of  the  initial  contact  surface.  In  order  to  make  this  testing 
time  as  long  as  possible  it  is  necessary  to  ensure  that  this  wave  system 
is  swept  through  the  test  section  in  as  short  a  time  as  possible^®)  The 
governing  factor  is  the  strength  of  Sg  which  for  ideal  starting  must 
be  a  sound  wave.  By  placing  a  diaphragm  at  the  nozzle  entrance  the 
test  section  can  be  evacuated  to  a  pressure  much  below  the  pressure 
in  the  uniform  tube  which  will  in  turn  strengthen  the  initial  shock 
wave  when  it  enters  the  nozzle.  Therefore  the  velocity  difference  across 
Sg  will  be  reduced  and  in  the  limit  Sg  will  become  a  sound  wave.  If 
the  initial  pressure  in  the  nozzle  is  decreased  still  further  expansion 
waves  will  be  formed  moving  towards  the  nozzle  entrance.  However  since 
expansion  waves  move  at  the  local  speed  of  sound  relative  to  the  flow 
there  is  no  gain  in  reducing  the  pi  ssure  ratio  across  the  nozzle 
diaphragm  below  that  value  required  to  make  Sg  a  sound  wave.  Fig.  18 
gives  the  pressure  ratio  required  across  the  nozzle  diaphragm  to  ensure 
y  r'V.-et  starling  for  a  range  of  test  section  Mach  numbers.  This  figure 
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indicates  that  an  extremely  large  pressure  ratio  is  required  to  obtain 
a  perfect  start  for  Mach  numbers  greater  than  10, 

Figure  19  shows  the  effect  of  reducing  the  pressure  ratio  across 
the  nozzle  diaphragm  thus  giving  an  imperfect  start.  In  this  case  the 
wave  S2  (Fig.  16)  will  stand  in  the  nozzle  and  suggests  the  possibility 
of  placing  a  plenum  chamber  downstream  of  the  test  section.  This 
chamber  will  create  an  expansion  wave  which  will  then  move  upstream, 
weaken  the  shock  Sg  and  pull  it  through  the  test  section.  A  schematic 
drawing  of  the  arrangement  is  shown  in  Fig.  20  while  Fig.  21  indicates 
the  pressure  ratio  required  of  the  plenum  chamber  for  a  range  of  hyper¬ 
sonic  Mach  numbers  in  the  test  section. 

A  plenum  chamber  may  also  be  employed  when  the  shock  Sg  is  not 
standing  in  the  test  section.  In  this  case  the  expansion  waves  from 
the  plenum  will  further  accelerate  the  starting  process  and  thus  in¬ 
crease  the  test  time.  If  the  pressure  in  the  nozzle  and  test  section 
were  low  enough  to  ensure  a  perfect  start,  addition  of  a  plenum  chamber 
would  not  improve  starting,  however,  the  addition  of  a  plenum  is  advan¬ 
tageous  in  another  way.  For  a  perfect  start,  S-j_  (Fig.  16)  passes  through 
the  test  section  at  a  high  velocity  and  a  long  test  section  is  required 
to  prevent  the  reflected  wave  returning  and  destroying  the  flow  before 
the  arrival  of  the  original  contact  surface,  C^.  By  adding  a  plenum 
chamber  at  the  end  of  the  test  section  this  reflected  shock  can  be 
greatly  weakened.  Therefore,  by  using  a  plenum  of  large  cross-sectional 
area  and  shock  length  the  facility  may  be  made  more  compact  and  even 
reduce  the  pumping  capacity  required. 
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6.3  THE  REFLECT?!)  METHOD 

If,  instead  of  the  straight  expansion  nozzle,  the  driven  section 
is  terminated  by  a  regular  convergent-divergent  hypersonic  nozzle  a 
reflected  type  shock  tunnel  results.  At  the  hypersonic  Mach  numbers 
of  interest  the  area  ratio  is  so  large  and  the  throat  area  so  small 
that  the  nozzle  will  act  as  a  solid  plate  end  the  initial  shock  wave 
will  be  completely  reflected.  This  means  that  the  gas  following  the 
shock  will  be  compressed  a  second  time  and  brought  to  rest  creating 
a  high  pressure,  high  temperature,  stagnant  reservoir  of  air  to  be 
expanded  through  the  nozzle.  The  reflected  shock  governs  the  duration 
of  steady  flow,  interacting  with  the  contact  surface  and  returning  to 
the  nozzle.  A  wave  diagram  and  sketch  of  this  method  is  shown  in  Fig. 

22. 

The  reflected  method  has  the  advantage  that  hypersonic  flow  is 
initiated  from  a  reservoir  with  almost  constant  supply  conditions. 
Therefore,  the  problems  of  shock  attenuation,  which  would  constantly 
vary  the  test  conditions,  and  boundary  layer  build  up  behind  the 
initial  shock  wave  are  not  present.  For  the  same  temperature  require¬ 
ments  the  initial  shock  may  be  weaker  than  in  the  nonreflected  case, 
or  conversely,  higher  stagnation  temperatures  are  obtained  for  the  same 
shock  strength.  The  testing  time  for  the  same  Mach  number  is  also 
slightly  increased  due  to  the  decrease  in  required  initial  shock  strength 
The  starting  problems  with  this  type  of  tunnel  are  essentially  the 
-is*  that  of  the  nonreflected  type  and  a  diaphragm  at  the  nozzle  is 
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desirable  in  order  to  .reduce*  the  starting  time.  In  the  nonreflected 
case  the  entrance  Mach  number  jnay  be  as  high  as  1.89,  based  on  perfect 
fluid  theory,  however,  for  the  reflected  type  the  Mach  number  will  be 
unity  at  the  throat.  Fig.  23  presents  the  pressure  ratio  required  across 
the  nozzle  diaphragm  over  a  range  of  test  section  Mach  numbers.  Comparing 
this  figure  with  Fig.  19  it  is  seen  that  there  is  a  reduction  in  the 
pressure  ratio  required  for  starting  in  the  reflected  case. 

6.4  THE  TAILORED-INTERFACE  MODIFICATION 

The  test  time  in  the  reflected  shock  type  of  tunnel  is  limited  to 
the  time  interval  between  the  arrival  of  the  initial  shock  wave  at  the 
nozzle  and  the  return  of  the  secondary  wave  generated  by  the  shock- 
contact  surface  interaction.  A  considerable  increase  in  testing  time 
is  possible  if  the  reflected  shock  can  be  made  to  pass  through -the 
contact  front  without  producing  any  additional  waves.  This  type  of 
tunnel  derives  its  name  from  this  impedance  matching  or  "tailoring", 
and  two  examples  are  illustrated  in  Fig.  24.  If  a  second' convergent- 
livergent  nozzle  is  placed  at  the  initial  diaphragm  station  then  the 

supersonic  portion  of  the  rarefaction  wave  can  be  eliminated  and  the 

» 

reflected  shock  will  traverse  the  full  length  of  the  driven  section 
before  returning  to  the  downstream  nozzle.  In  this  case  the  driver 
section  must  be  long  enough  to  ensure  that  the  reflected  rarefaction 
wave  0oc3  not  reach  the  downstream  nozzle  before  the  reflected  shock 
wave  returns.  One  disadvantage  of  the  method  is  that  the  area  ratio 
of  the  nozzle  at  the  initial  diaphragm  station  must  be  changed  for 
each  Mach  number.  & 
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If  the  reflected  shock  wave  passeG  through  the  contact  front 
without  creating  any  additional  waves  the  driver  and  driven  gases 
must  be  at  equal  pressure  and  velocity  behind  the  reflected  shock. 

As  the  driver  and  the  driven  gases  already  are  at  equal  pressure  and 
velocity  ahead  of  the  reflected  shock  wave,  the  pressure  ratio  and 
velocity  difference  across  the  reflected  shock  are  the  same  for  the 
two  gases.  This  may  be  stated  mathematically  with  the  aid  of  Fig.  24. 
The  conditions  for  nonreflected  waves  are  as  follows.^ 

£1- P 

*  p. 

ax  a , 

In  the  second  of  these  two  equations  we  have  a  Mach  number  and 
the  speed  of  sound.  In  this  case  the  strength  of  the  incident  and 
transmitted  shocks  are  equal  only  if  expressed  in  terms  of  their 
pressure  ratios.  The  change  in  velocity  across  the  shock  may  be 
expressed  in  the  following  manner. 


Combining  (6:4;l)  and  (6:4: 2)  gives 

.  V,  {*,-')(  IjJjz IP) 

ol  a ; 


(6:4:1) 


(6:4:2) 


(6:4:3) 
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This  equation  must  "be  satisfied  whenever  a  shock  passes  through  a  contact 
surface  without  creating  a  reflected  wave.  For  equal  values  of  V  the 
equations  then  depend  solely  on  the  speeds  of  sound  in  states  (2)  and 
(3).  This  would  be  the  case  for  a  i^/Air  combination.  Therefore  knowing 
the  initial  pressure  p]_,  the  shock  vave  Mach  number  and  the  driver  gas 
then  the  required  initial  pressure  and  temperature  of  the  driver  gas 
can  be  calculated. 

Figures  25  and  2 6  give  the  initial  diaphragm  pressure  ratio  and 
the  initial  driver  gas  temperature  required  for  tailoring  over  a  range 
of  test  Mach  numbers  for  helium  and  hydrogen  drivers.  It  is  seen  from 
Fig.  26  that  helium  at  room  temperature  will  tailor  at  a  Mach  number 
of  7.5,  hydrogen  at  10.6.  In  order  to  use  the  tailored  interface  tech¬ 
nique  at  higher  Mach  numbers  the  driver  gas  must  be  heated.  Figure  27 
gives  the  variation  of  flight  Mach  number  with  shock  Mach  number  for 
these  conditions.  Another  disadvantage  of  the  steady  configuration, 
indicated  by  the  figures,  is  the  increased  pressure  and  temperature  of 
the  driver  gas  required  for  this  modification.  However  the  steady 
tailored  interface  method  produces  test  times  three  times  greater  than 
the  unsteady  modification  and  twenty-five  times  the  nonreflected  method. 

A  comparison  of  testing  times  is  shown  in  Fig.  28. 

6.5  REAL  GAS  EFFECTS  IN  A  HYPERSONIC  SHOCK  TUMEL 

In  Section  5  the  real  gas  effects  in  the  production  of  strong 
shocks  were  considered.  A  brief  consideration  will  be  given  here  of  the 
real  gas  effects  on  the  design  of  a  hypersonic  shock  tunnel,  using  on 
expansion  nozzle.  Figure  29  indicates  the  variation  of  test  section 
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Mil'll  numbers  for  a  r .-mg'-  shock  Mach  numbers  .for  real  ,  calculated 

for  an  equilibrium  tost  section  s tat lo.  temperature  of  2l8°K,  corresponding 
to  tlie  isothermal  atmosphere .  The  curve  for  the  ideal,  gas  case  is  also 
shown  for  comparison,  and  is  approximately  C  per  cent  lower  than  the  real 
gas  results.  The  density  variation  is  shown  in  Fig.  30  and  the  real  gac 
values  are  decreased  compared  to  the  ideal,  gas.  Figure  31  presents  the 
area  ratio  variation  and,  as  would  be  expected  from  the  density  results, 
the  real  gas  area  ratio  is  considerably  increased.  In  the  reflected  type 
the  temperatures  ahead  of  the  nozzle  may  be  high  enough  for  dissociation 
to  occur.  When  the  gas  is  expand'1,  through  the  nozzle  the  cooling  rate 
vrill  be  very  high  and  some  of  the  constituents  such  as  NO  may  be  frozen. 
Another  problem  encountered  using  the  reflected  method  is  the  loss  of 
energy  by  radiation  from  the  hot  gas  behind  the  reflected  shock,  vhich 
becomes  more  important  as  testing  times  arc  increased.  These  increased 
test  times  are  possible  using  the  -tailored  interface  modification  and 
so  the  radiation  losses  will  have  to  be  considered  in  a  shock  tunnel 
of  this  type. 

0.6  REYNOLDS  NUMBER  AND  STAGNATION  TEMFER.'VfURE  SI.MNLATION 

f 

In  order  to  do  hypersonic  testing  it  is  desirable  to  simulate  free 
flight  conditions  as  closely  as  possible.  This  means  simulating  Mach 
number,  stagnation  temperature  and  Reynolds  number.  In  the  straight 
shock  tube  it  ic  possible  (.0  simulate  the  temperature  and  Reynolds  number 
but  the  Kno/i  number  is  limited  to  a  value  of  three  as  ohovn  in  Fig.  13. 


J’i.e  US'!  '.nonunion  noiw.les  •«,  the  downs  rn-or.-n  end  of’  the  driven  section. 


f 
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enables  the  flow  to  expand  to  a  high  Mach  number  but  in  doing,  so  may 
limit  the  other  two  conditions.  The  correct  Reynolds  number  and  stagnation 
temperature  are  controlled  by  the  initial  condition  before  the  test  ga3 
enters  the  nozzle  and  thereby  the  strength  of  the  initial  shock  wave. 
Figures  32  and  33  show  the  Reynolds  number  variation  with  test  Mach 
number  for  u  perfect  gasf(  H  -  1.4). and  a  test  section  static  temperature 
of  518. 4°R  for  a  range  of  pressures  behind  the  initial  shock  wave.  Due 
to  the  decrease  of  density  caused  by  the  real  gas  effects  (Fig.  30 )  the 
actual  Reynolds  number  will  be  lower  than  those  shown.  For  a  given 
3fcatic  pressure  in  the  driven  section,  higher  Reynolds  numbers  are 
realized  with  the  nonreflected  technique.'  Even  though  the  Reynolds 
numbers  shown  here  are  lower  than  those  possible  in  region  (2)  of  a 
straight  shock  tube  this  may  not  be  important  as  the  Reynolds  number 
of  a  body  traveling  at  a  hypersonic  Mach  number  (pr  20),  at  high  altitude 
(100,000  ft.),  is  relatively  low  (2  X  106/ft.).  Therefore  the  attain¬ 
ment  of  a  high  Reynolds  number  will  not  be  difficult  for  this  class  of 
problems . 


7  CONCLUDING  REMARKS 

This  report  has  attempted  to  outline  the  basic  theory  of  the  shock 
tube  end  describe  most  of  the  modifications  to  the  basic  tube  that  have 
boon  devised  for  hypersonic  testing.  Detailed  calculation  and  theory 
have  been  omitted  but  these  may  be  found  in  the  many  references  cited. 
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APPENDIX  A 

NOTATION 


Figure  1  serves  to  define  the  flow  regions  in  the  basic  shock 
tube.  The  notation  Air/Air,  He/Air  and  Hg/Air  denotes  a  gas  com¬ 
bination  across  the  diaphragm,  the  left-hand  symbol  representing 
the  gas  in  the  driver  section,  the  right-hand  symbol  the  gas  in  the 
driven  section. 


A 

Aid 

CU 

D 

Eij 

e 

L 

Ms 

M 

P 

*o 

p« 

Pl*d 


t 

T 


Velocity  of  sound 
Area 

Dimensionless  speed  of  sound  ratio  (i.e.  Aq4  =  a-^/a^) 

Dimensionless  rarefaction  wave  speed 

Diameter 

Dimensionless  energy  ratio  (CvT)j/  (CyT)^ 
Equivalence  factor  (Eq.  4;3;2) 

Length 

Shock  Wave  Mach  number  based  on  speed  of  sound  of  gas 
into  which  it  travels 

Mach  number 

Pressure 

Overall  pressure  ratio  (see  Section  4.2) 

Dimensionless  pressure  ratio  a...  *b/n> 

Pressure  in  region  (4),  Fig.  4  for  double  diaphragm 
shock  tube  (Sect.  4.2) 

Pressure  in  region  (4)  for  single  diaphragm  shock  tube 
(Sect.  4.2) 

Time 

Temperature 
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,  TU 


u 


Uij 

X 


* 


Dimensionless  temperature  ratio  (i.e.  Ta/Tj.) 

Temperature  in  region  (^i  )<  Fig.  *1  for  double  diaphragm 
shock  tube  (Sect,  h.2) 

Temperature  in  region  (it)>  for  single  diaphragm  shock 
tube  (Sect.  U.2) 

Particle  velocity 

Dimensionless  particle  velocity  (l.e.  Ug/aj.) 

Position  along  the  shock  tube  measured  from  the  diaphragm 
Ratio  of  specific  heats  (C  /cv  ) 

Jr 


Hi 

v-i 

hi 

2V 

Dimensionless  density  ratio  (i.e.  pg/  pj_) 

Kinematic  viscosity 
Density 

Contact  surface  moving  to  the  right 

Backward  facing  rarefaction  wave  (particles  enter  from  left) 
Forward  facing  shock  wave  (particles  enter  from  light) 


Quantities  not  included  here  are  defined  in  the  text  or  on  the 
figures  in  which  they  appear. 
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10 


‘21 


Figure  2  -  Variation  of  Shock  Pressure  Ratio,  with  Diaphragm  Pressure  Ratio, 


Pj^,  7  m  l»^-« 
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FIG.  4  -  The  Double  Diaphragm  Shock  Tube  -  Reflected 
Shock  Type 


FIG.  '}  -  The;  Double  Diaphragm  Shock  Tube  -  Unsteady 
Expansion  Type 
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HO.  7  -  Shock  Tube  with  a  Convergent  Trnnciton  Section.  The 

dnr.hed  linen  in  the  diaphragm  section  denote  a  convergent 
-divergent  geometry  with  minimum  area  at  3b* 
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-  Performance  Comparison  of  a  Variable  Geometry  and 
Simple  Shock  Tube,  7  *  l„k. 


^3 


Figure  8 


initial  pressure  ratio: suffer  to  air 

Figure  9  -  Performance  of  a  hydrogen-air  shock  tube  with  an  over-all  pressure  ratio  of  10 
showing  the  effects  of  varying  the  pressure  and  molecular  weight  of  the  buffer 
(Ref.-  12) 
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,ture  Batlo 


Shock  Mach  Number,  M# 

Figure  12  -  Density  Ratio  Variation  across  the  Initial 
8hock  Wave. 
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Shock  Mach  Number.  M 

B 

Figure  13  -  Variation  of  Flow  Mach  Number  boa 
Initial  Shock  Wave 
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Figure  15  -  Overall  area  ratio  vs.  flight  Mach  number  (7  •  1*4. 
no  reflection)  (Ref.  36) 
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17  -  The  Cornell  Aeronautical  Laboratory  11-  By  15-  Inch  Hypersonic  Shock  Tunnel 
(Ref.  12) 
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STEADY  MODIFICATION 


Figure  2k  -  Wave  diagrams  cf  shock  tube  modifications  for  Increased 
testing  time.  (Ref.  12) 
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FUORT  MACH  HUMOER 

Figure  25a  -  Required  di.apbre.gui  pressure  ratio  for  tailoring  at  various 
flight  Mach  numbers  (unsteady  configuration) .  (Ref,  39) 
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Figure  25b  -  Required  diaphragm  pressure  ratio  for  tailoring  at 
various  flight  Mach  numbers  (steady  configuration; 
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Figure  26b  -  Required  initial  driver  gas  temperature  for  tailoring  at 
various  flight  Mach  numbers  (steady  configuration). 

(Ref.  39) 
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SHOCK  MACN  NUMBEK 

Figure  27  -  Variation,  of  flight:  Mach  number  with  shock  .Mach 
numbers,  (steady  configuration) , (Ref .  39) 
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Figure  28a  -  Testing  time  for  four  modifications .  (Ref .  37 ) 
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FLIGHT  MACH  NUMBER 

Figure  28b  -  Nominal  testing  time  of  hypersonic  shock  tunnel  modifications 
based  on  equilibrium  real  gas  calculations,  (Ref.  12) 
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DRIVER  SHOCK  MACS  HUMBER,  H, 

Figure  29  -  Variation  cf  Test  Section  M.ach  number  with  driver 
shock  Mac: •  number.  (Ref.  14) 
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Figure  30  -  Density  in  test  section  of  hypersonic  shock  runnel  vs. 
test  section  Mach  number .  (Ref.  14) 
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TEST  SECTION  MACH  NUMBER,  Mf 


Figure  31  -  E".parjoioa  or  the  hypersonic  shocl;  tunnel 
section  Mach  number’  (ner'  Ik  ) 
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•'i.'-'it  NUMBER  Mf 

Figure  33  -  Perforinauce  of  the  reflected  method  T.  »  516,4°R, 

*  391  R,  R  ■  1716  FT^/SEC2  °R;  «  static  pressure 

behind  the  reflected  shock  (Lr./lN2)  (Ref.  37) 


quantities  for  very  strong  shock  vave6. 
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*0 

10 

102 

103 

101* 

105 

m 

1.64 

2.66 

4.39 

7.20 

11.79 

TABUS  2 

MAXIMUM  8H0CK  MACH  NUMBERS  FOR  MULTIPLE  DIAPHRAGM  SHOCK  TUBES  (Ref ,  9) 


Driver  -  Driven 

Ali/Aj^ 

Oaa  Combinations 

1 

19 

2.25 

00 

*21*2 

6.18  * 

6.34 

6.47 

6.76 

VA 

7.48 

.7.69 

7.84 

8.20 

*4i  *  1 

He/«2 

10.9 

11.3 

11.6 

12.5 

w* 

12.8 

13.3 

13.7 

14.7  . 

%/h2 

22.6 

.23.2 

23.6 

24.7 

|| 

27.5 

.28.3 

28.8 

30.1 

■1 

'  He/N2 

14.8 

15.4 

15.9 

17.1 

He/A 

18.0 

18.8 

19.3 

20.8 

mm 

Ha/Na 

31.9 

32.8 

33.4 

34.9 

H 

Hp/A 

38.9 

40.0 

40.7 

42.6 

TABLE  3 

ICVKDSM  VALUES  0?  SHOCK  MACH  HUMBER,  Ms  FOR  P;4l  -*•  00  (Ref .  10 ) 
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